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Seismology of Procyon A: determination of mode frequencies, 
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Abstract. The F5 IV-V star Procyon A (a— CMi) was observed in January 2001 by means of the high resolution spectrograph 
SARG operating with the TNG 3.5 m Italian telescope (Telescopio Nazionale Galileo) at Canary Islands, exploiting the iodine 
cell technique. T he time-series of abo ut 950 spectra carried out during 6 observation nights and a preliminary data analysis 
were presented in IClaudietal.ld2005h . These measurements showed a significant excess of power between 0.5 and 1.5 mHz, 
with =! lms" 1 peak amplitude. Here we present a more detailed analysis of the time-series, based on both radial velocity and 
line equivalent width analyses. From the power spectrum we found a typical p-mode frequency comb-like structure, identified 
with a good margin of certainty 11 frequencies in the interval 500 - 1400 /jHz of modes with I = 0, 1,2 and 7 < n < 22, and 
determined large and small frequency separations, Av = 55.90 ± 0.08 /iHz and tSvm = 7.1 ± 1.3yuHz, respectively. The mean 
amplitude per mode (Z = 0, 1) at peak power results to be 0.45 ± 0.07 ms~', twice larger than the solar one, and the mode 
lifetime 2 ± 0.4 d, that indicates a non-coherent, stochastic source of mode excitation. Line equivalent width measurements do 
not show a significant excess of power in the examined spectral region but allowed us to infer an upper limit to the granulation 
noise. 
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1. Introduction 

Procyon A (a - CMi, HR 2943, HD61421), hereinafter sim- 
ply called Procyon, is a F5 IV-V star with V = 0.363 at a 
distance of 3.53 pc in a 40 y period visual binary system, the 
companion being a white dwarf more than 10 mag fainter. By 
adopting the very accura te parallax measured b y HIPPARCOS, 
n = 285.93 + 0.88 mas. lAUende Prieto et al.l d2002l) derived a 
mass of 1 .42 + 0.06 M , a radius of 2.07 1 + 0.02 R , and a grav- 
ity logg = 3.96 + 0.02. 

Procyon, owing to its proximity and brightness, 
eady attracted attention of stellar seismolo- 



However, data f rom the Canadian MOST satellite 
(Matt hews et al. 2004) show no sig nificant power excess in the 
same spectral region. To this regard lBedding et al. I (12005b sug- 
gested that the most likely explanation for the null detection 
could be a do minating non-stellar no i se so urce in the MOST 



data, though Reg ulo & Roca Cortes I (120051) claimed for the 



has 



gists (Brown et al. 1991 ; Guenther & Demarque 


1993 


Barban et al. 1999; Chaboyer et al. 


1 19991 iMartic etal. 


1999 


2004; lEggenberger et al. 


2004), causing an intense debate 



presence of a signal in these data. 

In Paper I we presented high precision radial velocity mea- 
surements carried out during 6 observation nights by means 
of the high resolution spectrograph SARG operating with 
TNG 3.5 m Italian telescope (Telescopio Nazionale Galileo) at 
Canary Islands. The data showed an excess of power between 
0.5 — 1.5 mHz with a large separation of abou t 56 + 2pHz, 
in agr e ement with p r evious measurements of iMosser et al. I 



among the scientists. T he excess of power in the range 



(1998), Barba n et al 



0.5 - 1 .5 mHz found bv iBouchy et al. I d2004. IMartic et al. 



(2004), Egge nbergeret al.l (12004 . and IClaudi et al 



1 (fl99 
(1200- 



Marti cetal 



i of 5 
(11999 



who found 



19991 120041) and 
values in the range 



(2005) 



hereinafter referred as Paper I, seems to have a stellar origin 
and is consistent with a p-mode comb-like pattern. 



Send offprint requests to: S. Leccia; e-mail: leccia@na.astro.it 



Eggenberger et al 
53 - 56 //Hz. Here we deal with an improvement of the pre- 
liminary analysis presented in Paper I concerning an accurate 
time-series analysis of both radial velocities and line equiva- 
lent width measurements, from which mode frequencies, am- 
plitudes and lifetimes are determined and constraints on stellar 
granulation noise are deduced. 
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2. Observations 

SARG is a high reso lution cross dispersed echelle spectrograph 
dGratton et al. 11200 lb which operates in both single object and 
long slit (up to 26") observing modes and covers a spectral 
wavelength range from 370 nm up to about lOOOnm, with a 
resolution ranging from R - 29, 000 up to R = 164, 000. Our 
spectra were obtained at R — 144, 000 in the wavelength range 
between 462-792 nm. The calibration iodine cell works only in 
the blue part of the spectrum (462 - 620 nm) that has been used 
for measuring Doppler shifts. The red part of the echelle spec- 
trum of Procyon has been used for measuring line equivalent 
widths of absorption lines sensitive to temperature. During the 
observing period we collected about 950 high signal-to-noise 
ratio (S/N) spectra with a mean exposure time of about 10 s 
(see Paper I for more details). The red part spectrum has been 
calibrated in wavelength by using a Th-Ar lamp. The analysis 
of both the blue and red parts of the spectrum was performed 
by using the IRAF package facilities. 



3. Radial velocity measurements and analysis 

Radial velocities have been dete rmined by means of the 
AUSTRAL code fendl et al. Il200oh which models instrumen- 
tal profile, stellar and iodine cell spectra in order to measure 
Doppler shifts (see also Paper I). This code also provides an 
estimate of the uncertainty in the velocity measurements, cr,. 
These values have been derived from the scatter of velocities 
measured from many (^ 100), small (^ 2 A) segments of the 
echelle spectrum. In order to perform a weighted Fourier anal- 
ysis of the data, we firstly verified that these cr,- values reflected 
the noise propertie s of the velocity measurements, following 
Butl er et al. ( 2004) approach. The high frequency noise in the 
power spectrum (PS) well beyond the stellar signal, reflects the 
properties of the noise in the data, and since we do expect that 
the oscillation signal is the dominant cause of variations in the 
velocity time series, we need to remove it before to analyze 
the noise. We do that iteratively by finding the strongest peak 
in the PS of the velocity time-series and subtracting the corre- 
sponding sinusoid from the time-series. This procedure is car- 
ried out for the strongest peaks in the oscillation spectrum in 
the frequency range 0-2 mHz, until the spectral leakage into 
high frequencies from the remaning power is negligible. Thus 
we have a time series of residual velocities, r, , that reflects the 
noise properties of the measurements. We then analyze the ra- 
tio r,/crj, which is expected to be Gaussian-distributed, so that 
the outliers correspond to the suspected data points. The cu- 
mulative histogram of |r,7cr,j is shown in the upper panel of 
Fig- HI The solid curve shows the cumulative histogram for the 
best-fit Gaussian distribution. A significant excess of outliers 
is evident for |r,/cr,| > 2. The lower panel of Fig. Q] shows the 
ratio of the values of the observed points to the corresponding 
ones of the Gaussian curve, i.e. the fraction / of data points 
that could be considered as "good" observations, namely those 
which are close to the unity. The quantities w, = 1 / (ajf) have 
been adopted as weights in the computation of the weighted PS 
shown in Fig.|2l where the most prominent peak has an ampli- 
tude of ^ lms 4 . 




r/a 



Fig. 1. Upper panel: cumulative histograms of |r,/cr,| for SARG 
data. The diamonds show the observed data, and the solid curve 
shows the result expected for a Gaussian-distributed noise. 
Lower panel: ratio of the observed to the expected histograms, 
indicating the fraction of "good" data points. 

3.1. Search for a comb-like pattern 

In solar-like stars, p-mode oscillations of low harmonic de- 
gree, I, and high radial order, n, are expected to produce in the 
PS a characteristic comb-like structure with mode frequencies 
v„j reasonably well app roximated by the Tassoul's simplified 
asymptotic relationship (iTassoul l ll980h : 



Av(n + 1/2 + e) - 1(1 + l)<5v 02 /6 



(1) 



where Av =< v„j - v„-ij > and <5v 02 -< >7?,o - v„_i >2 > are the 
average large frequency separation and small frequency separa- 
tion for / = 0, 2, respectively, and e is a constant of the order of 
unity, sensitive to the sound speed near the surface layers of the 
star. Av reflects the average stellar density, while 5vq2 is largely 
determined by the conditions in the core of the star and reflects 
its evolutionary state. In order to estimate the large frequency 
separation Av in the region of the excess of power detected in 
the PS, we applied the comb-response method, that is a gener- 
alization of the PS of a PS and consequently allows us to search 
for any regularity in a spectral pattern. The comb-re sponse 



function, CR(Av), is defined by Kjel dsen et al. I (119951) as the 



following products of PS calculated for a series of trial values 
Av: 

CR(Av) = PS(f_i /2 )PS(f +I / 2 )PS(?-i)PS(f + i) x 
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Fig. 2. The power spectrum of the weighted data. An excess 
of power around 1 mHz is evident. The inset shows the power 
spectrum of the window function for a sine-wave signal of am- 
plitude 1ms -1 , sampled in the same way as the observations. 
The power units are the same as those of the main figure. 



600 800 1000 

Frequency (,".Hz) 



Fig. 3. The first-order spacing calculated over several fre- 
quency ranges, as determined from the comb-response analy- 
sis. The line shows the mean value of 55.7 //Hz, with a standard 
error of + 1.4 //Hz. 



pS(f_ 3/ 2)PS(f + 3/2)PS(f-2)PS(? + 2)] 



1/2 



(2) 



where = v max +£Av, with v max the frequency of the strongest 
peak in the PS. A peak in the CR at a particular spacing Av 
indicates the presence of a regular series of peaks in the PS, 
centred at the central frequency v max - tentatively assumed to 
be an oscillation mode - and having a spacing Ay/2. In or- 
der to select the central frequencies v's for the CR analysis 
we evaluated the white noise in the PS between 2.5 - 3.0 mHz 
and adopted as v max 's the frequencies of the peaks in the PS 
between 0.5 - 1.5 mHz with S/N > 4. For each central fre- 
quency v m ax we searched for the maximum CR in the range 
20 < Av < 80 //Hz. Figure [3] reports the Av's determined from 
the CR analysis, that gives a mean value Av = 55.7 + 1.4 //Hz, 
where the error is the standard deviation, and Fig.|4]the cumula- 
tive CR function computed as the sum of all the individual CR 
functions referring to each central frequency v max , that gives 
Av = 56 + 1 //Hz, where the error is the FWHM of the peak. 



Fig. 4. The cumulative comb-response obtained as the sum of 
the individual comb-responses for each central frequency vo. 
The peak is centred at Av = 56 + 1 /(Hz. The error is the FWHM 
of the peak. 



Frequency MOD 56 (>Hz) 



Fig. 5. The echelle diagram as derived from the "standard" ex- 
traction of frequencies. 



3.2. Oscillation frequencies and mode amplitude 

By assuming the 56 /(Hz regularity as genuine, we attempted 
to identify the oscillation frequencies directly from the PS. In 
order to extract the frequencies we used the following proce- 
dure ("standard" extraction): i) we found the largest peak, in 
the frequency range 0.5 - 1 .5 mHz, then subtracted the corre- 
sponding sine-wave function from the time-series, finally re- 
computed the PS; ii) we repeated such a procedure for the new 
largest peak in the same range of frequency, the procedure be- 
ing stopped when there were no peaks with amplitude larger 
than 0.29 m 2 s -2 , namely larger than 3<x, where <x is the white 
noise evaluated in the PS between 2.5 - 3 mHz. The extracted 
frequencies are plotted in the echelle diagram shown in Fig. [5] 
Owing to single site observations, we need to take into account 
the daily alias of 1 1 .57 //Hz (see the spectral window in Fig. |3J 
in order to identify the modes in the echelle diagram; in fact, we 
cannot know a priori whether or not the frequencies, selected 
by the iterative procedure described above, are the correct ones 
or just aliases. Therefore, instead of shifting the frequencies 
by the daily alias, that, for each shift, might lead to an arbi- 
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Table 1. Prominent peaks in the power spectrum of Procyon 
identified as mode oscillation frequencies (in //Hz) in terms of 
n and I. The values of n are deduced from the fit to the asymp- 
totic relationship ([TJ. The frequency indicated in parenthesis 
was excluded from the fit because of its large deviation from 
the asymtotic relationship, though the related peak in the power 
spectrum is a high S/N one. 
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Fig. 6. The echelle diagram as derived from the "modified" ex- 
traction of frequencies. 



trary identification of the modes, we prefer to use a procedure 
("modified" extraction) which relies on two assumptions only: 
i) the largest peak in the range of interest is a true pulsation 
mode; ii) the large separation value is 56 //Hz. We operated 
in the same way as in the point i) of the "standard" extrac- 
tion, then defined a new searching spectral region centred at a 
distance of 56 //Hz from the first component, the extension of 
which was 2 //Hz, a suitable frequency span related to spectral 
resolution. If this region contained a peak with an amplitude 
larger than 3<x, we identified this peak as the second compo- 
nent to be cleaned, and recomputed the PS. This process was 
repeated, shifting frequency back and forward by multiples of 
56 //Hz, and once the whole region of interest had been cov- 
ered, the remaining largest peak was identified and the process 
restarted with a new set of searching regions separated again by 
56 //Hz, the procedure being stopped when all remaining peaks 
were bel ow the 3<x threshold. This is the same procedure as that 
used by Kjeld sen et al. I dl995 ) for 77 Boo, that allows to deter- 
mine the mode frequencies except for the daily alias that could 
cause a shift of + 11. 57 //Hz. The frequencies extracted with 
this procedure are shown in Table [TJ and the related echelle di- 
agram in Fig. [6] The large and small frequency separations, 




Frequency MOD 55.9 fiHz 



Fig. 7. The mean power folded at the large separation frequency 
as a function of frequency modulo of 55.90 //Hz. The positions 
of the modes with / = 0, 1 , 2 are indicated together with the 
positions of the daily side-bands. 



and the constant e have been determined by means of a least 
square best fit of frequencies with the asymptotic relationship 
O, obtaining: Av = 55.90 + 0.08 //Hz, 6v 02 = 7.1 + 1.3 //Hz, 
e — 1.913 + 0.025, with a rms scatter per mode of 1.23 jjHz. In 
Table [JJ we also give a possible mode with I = 1 and n — 12 
at 808.6 /iHz, however this mode was excluded from the calcu- 
lation of the fit because of its large deviation from the asymp- 
totic relation, though the related peak in the power spectrum is 
a high S/N one. Once determined Av, it is possible to construct 
a folded spectrum where the power in the PS is folded at the 
large separation, by using the frequency modulo of 55. 90 //Hz. 
The result is shown in Fig. |7]where the mean power as a func- 
tion of frequency modulo is plotted, and the positions of modes 
with I — 0, 1 , 2 are given together with the positions of the daily 
side-bands, from which the detection of a p-mode structure in 
Procyon is evident. 

The amplitude of individual modes can be estimated from 
the power concentrated in the PS, on assuming that, owing to 
the full disc observations, only I = 0, 1 , 2, 3 modes are detected. 
It is then possible to estimate the amplitude per mode necessary 
for producing th e observed power leve l. Following the proce- 
dure described in lKieldsen et al. I (12005b we found a mean am- 
plitude of the p-mode peak power, for the modes with / = 0, 1 
in the frequency interval 0.6- 1 .4 mHz, of 0.45+0.07 m s _1 , that 
is a value twice larger than the solar one a nd in very goo d agree- 
ment with the amplitudes estimated by iBrown et al. I dl99ll) . 
Our velocity ampli tude can be transformed in intensity ampli- 
tude by using the lKjeldsen & Bedding I (119951) method, and ob- 
taining, at 5500 A, the value of 7.3 + 1.1 ppm per mode ( / = 0, 1) 
in good agreeme nt with the WIRE (|Bruntt et al. 112005 ) and up- 
per limit MOST ([Matthews et al. 120041) measurements. The re- 
sults are summarized in Fig. [8] where the velocity amplitude per 
mode is shown as a function of frequency together with the in- 
tensity scale and the equispaced positions of / = modes as 
derived from the relationship ([JJ in which the determined val- 
ues of Av and e were inserted. In Fig.[8]a dip at 1 mHz is clearly 
visible, that is a confirmation of a previously reported result of 
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Fig. 8. Velocity amplitude per mode (/ = 0, 1) as a function 
of frequency. The horizontal line at 0.45 ms~' represents the 
mean amplitude of the p-mode peak power in the frequency in- 
terval 0.6- 1.4 mHz. The scale of intensity is also given together 
with the positions of I = modes, as derived from asymptotic 
relationship (Q]i. 



Bedd ing & Kieldsen (2006), and whose presence is probabl y 
related to the damping rate of the modes (IHoudek et al. II 1999b . 

3.3. Comparison with other frequency determinations 

Our results are direc tly comparable with those of 
Eggenberger et al. ] d2004 . as given in their Table 3, where 
the identified frequencies are tabulated. We may also make 
a comparison with the r aw extracted frequencies given by 
Eggenberger et al. (2004) in their Table 2. We are not able 



to mak e a direct comparison with the results by Martic et al. 
(2004), since they only provide a list of identified mode 



frequencies that depend on the fitting to theoretical models. 



In order to perform the comparison with Eggen berger et al 
(2004) we apply the same leas t square best fit we appl ied to 
our frequency data to those of Eggenberger et al. I ([2004) and 
derive the parameters of the asymptotic relationship which re- 
sult to be: 6vq2 = 4.9 //Hz, s = 1.537. In this procedure 
the value of Av has been taken fixed to the value given by 
Eggenberger et al. (2004), Av = 55.5 //Hz. If one calculate the 
frequencies corresponding to the two asymptotic solutions, the 
one in the present paper a nd the one we obtain by fitting to the 
Eggenbe rger et al.l d2004l) data, we find that most of the calcu- 
lated frequencies, for I = 0, 1, agree within 3 //Hz. We find that 
I — 0, n = 9-24 modes from the Eggenberger et al. 



asymptotic relation correspond to / = 1, n = 8 -23 modes from 
the asymptotic relati on in the present study. W e also find that 
I = 1, n = 7 - 22 ( Eggenberger et al. d2004l) ) modes corre- 
spond to / = 0, n = 7-22 modes identifie d by us. We con- 



clude that we and Eggenberge r et al. I (120041) are detecting sig- 



nal from the same underlying p-m odes. If we look at all the 



Table 2. Comparison between the Eggenberger et al. I (|2004) 
extracted raw frequencies (in //Hz) and the asymptotic relation 
from the present study (Av = 55. 90 //Hz, <5vo2 =7.1 //Hz, s - 
1.913). 
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It is interesting to note that 8, out of 1 1, frequencies deter- 
mined by us match well, within a fe w //Hz, those listed in the 
Table 2 of Egg enberger et al. (2004), but with a different mode 
identification. Therefore we dete c ted si gnals from the same 
frequencies as Eggenberg er et al. 1 d2004 . but extracted the p- 
mode structure in a slightly different way, reaching two dif- 
ferent solutions, and therefore identifications, for the extracted 
frequencies. 

If we use our solution, Av = 55. 90 //Hz, <5vo2 = 
7.1 juHz, g = 1.913, to v erify to which extent it matches the 
Eggenbe rger et al. (2004) frequencies we find that 16 frequen- 
cies match our solution, 5 need a shift of +11. 6 //Hz, 4 need 
a shift of -11. 6 //Hz, an d only 2 do not ma t ch the solution. 



Otherwise, if we use the Eggenberger et al. I ([2004) solution, 



Av = 55.50//Hz, 6v 02 = 4.9 //Hz, s = 1.537, we find that 12 
frequencies match this solution, 8 need a shift of + 1 1 .6 //Hz, 
3 need a shift of - 1 1 .6 //Hz, and 4 do not match the solution . 



raw extracted frequencies given by Eggen berger et al. I (120041) This m eans that the frequencies detected bv lEggenberger et al 



in their Table 2, it appears of striking importance that most of 
those frequencies may be identified by using the asymptotic 
solution deduced from the presen t study. In table |2] we com- 
pare the Eggenber ger et al. I ( 20041) raw frequencies with the 
asymptotic relation from the present study (Av = 55.90 //Hz, 
(5v 02 = 7.1 //Hz, s= 1.913). 



( 2004) provide a slightly better fit to our solution than to their 
solution. 

This statement is supported by the fact that if we only 
look at those 12 frequ encies that fit, without shift, the 
Eggenberger et al. 1 d2004l) asymptotic rela tionship and com- 



pare them with the Eggenberg er et al. I ([2004) 16 raw frequen 
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cies that fit the solution found in the present study, we see that 
the scatter for the 16 frequencies is abou t 10% lower when 
we us e our asymptotic relation. Moreover, Egg enberger et al. 
(2004) identify the 50% of those 12 frequencies that fit its rela- 
tion as / = 2 modes, while the / = 1 modes should show instead 
the highest amplitude. 

However, since the two solutions provide basically the 
same frequencies, one should of course be cautious in claim- 
ing that one solution is far better than the other one. 



3.4. Mode lifetimes 

Solar-like p-mode oscillations in Procyon are presumably ex- 
cited stochastically b y the action of the surface convection 



( Houdek et al. I [19991) as in the case of the Sun's oscillations. 



Therefore it is interesting to give an estimate of the lifetime of 
modes for verifying the excitation mechanism that in turn de- 
termines their observed spectral width and amplitude. A coher- 
ent, over-stable, long-lived pulsation mode would concentrate 
all its power in a small frequency bin, whilst a stochastically, 
short-lived excited mode would spread its power over a larger 
frequency interval. In the case of the Sun the mode lifetime, 
at the peak amplitude of 3 mHz, is about 3 + 1 d, as estimated 
by fitting a model to the pow er di stribution for each m ode in 
the PS Chaplin etal.lll997h . The lChaplin et al. I d 1997b tech- 
nique requires observations for a significantly longer time than 
the mode lifetime, that is not necessarily consistent with our 
observing period of 6 nights. 

The method we use in the present a nalysis is based on the 



frequency scatter per identified mode (Bed ding et al. 1 12004 



Kjel dsen et al. 1 120051) . as calculated in Section 13.21 We ran 



8400 individual simulations in order to determine the rms scat- 
ter of the detected frequencies as a functio n of the mode life- 
time, by using the iDe Ridder et al. I d2006l) simulator, and ne- 
glecting rotational splittings both because Procyon is a slow ro- 
tator and the majority of modes used are radial. The simulations 
contain the mode correct amplitude, noise level, and number of 
detected modes out of the total number of modes contained in 
the p-mode spectrum, as determined from the asymptotic re- 
lationship (Q} in the observed frequency interval. The fit rep- 
resents a model that contains two components: the frequency 
spread, that is essentially the width of the Lorentzian profile, 
and the scatter of coherent oscillations, that corresponds to an 
infinite mode lifetime. The results are reported in Fig. [9] where 
the rms scatter of the 10 considered frequencies, out of the 1 1 
identified (see Section [3~2t . is plotted vs. mode lifetime. Since 
the mode frequency scatter determined in Section l3~2l results to 
be 1.23 //Hz, it is easily seen from Fig.[9]that the correspond- 
ing lifetime is 2.0 ± 0.4 d. The observed frequency scatter in 
Procyon is sensibly larger than that consistent with coherent os- 
cillations, ruling out this possibility, and the related mode life- 
times appear to be slightly shorter than those observed in the 
Sun, where the oscillations are excited by the surface turbulent 
convection. The present data strongly indicate that oscillations 
in Procyon are not coherent, but damped and re-excited by a 
mechanism similar to that operating in the Sun, as expected. 



10.0 




.0 10.0 

Mode lifetime (d) 



100.0 



Fig. 9. The rms scatter of the 10 identified frequencies as a 
function of mode lifetime obtained by several thousand sim- 
ulations. The intersection of the horizontal line of the observed 
scatter with the solid curve obtained by simulations indicates 
that the mode lifetime is 1 •97^91? d. The horizontal line of 
coherent oscillation scatter, corresponding to an infinite mode 
lifetime, is also indicated. The frequency spread reflects essen- 
tially the width of the Lorentzian profile. 



The results of the previous method for determining the 
mode lifetimes are corroborated by a different analysis we per- 
formed, based on the amplitude of the highest peaks in the PS, 
with the idea that also the amplitude of modes are affected by 
their lifetimes, in the sense that when the peak amplitude is 
smaller the mode lifetime is shorter, whilst when it is larger the 
mode lifetime is longer. However, since there is not a simple 
relationship between amplitude and lifetime, we used a data 
simulation, based on several thousand runs, that allows a direct 
comparison between the observed properties of the time-series 
and the properties of the oscillation modes in order to estab- 
lish a calibrate d relationship between the peak height and the 
mode lifetime (IDe Ridder et al. 1 12006). With respect to the pre- 
vious one, this method has the advantage of being independent 
of mode identification but the disadvantage of being less ac- 
curate. The result is reported in Fig. [TUl where we show the 
amplitude, with 1 cr error bars, of the 6 largest amplitude de- 
tected modes with / = 0, 1 , 2 as a function of the mode lifetime. 
In Fig.[l0]the intersection of the horizontal line corresponding 
to that of the observed mode of largest amplitude (see Fig. 13 
with the simulated data gives an average mode lifetime =* 2 d, 
in good agreement with the frequency scatter method but with 
a much larger uncertainty. It is interesting to note the expected 
trend of mode amplitudes as a function of their lifetimes. 

The effects of frequency scatter and mode amplitude on 
mode lifetime is evident in Figs. [TT] and [12J where simula- 
tions concerning 10 d continuous observation campaigns are 
reported for mode lifetimes of 2 d and 0.5 d respectively. In the 
above quoted two figures the upper panel shows the complete 
PS, while the lower one an enlarged portion of it from 0.8 mHz 
to 1 .0 mHz. From our simulations it is clearly seen that the fre- 
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0.80 _ 
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Mode lifetime (d) 



Fig. 10. The 

modes, for I 
bars indicate 



amplitude of the 6 largest amplitude detected 
= 0, 1 , 2 as a function of the mode lifetime. The 
the 1 <x error. 
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Fig. 12. The same simulation as that of Fig.[TT] but for a mode 
lifetime of 0.5 d. Note that the frequency spreads are larger and 
peak amplitudes are smaller than those of the modes seen in 
Fig. EH 

Table 3. Spectral line wavelengths used for EW measurements. 



Wavelength (A) 


Spectral lines 


6393.4 


Fel 


6436.7 


Cal 


6462.4 


Ca I / Fe I 


6546.0 


Fel 


6574.9 


Fel 


6592.7 


Fel 


6643.4 


Nil 


6677.8 


Fel 


6717.5 


Cal 



Fig. 11. Simul ated Procyon p-mode p ower spectrum obtained 
by using the iDe Ridder et aL (2006) simulator not including 
noise and rotational splittings. The length of the time-series is 
10 d and the mode lifetime is 2 d. 



quencies of long-living modes are much better defined, or less 
scattered, and have larger amplitudes than those of short-living 
modes. 



4. Line equivalent width measurements and 
analysis 

As already mentioned, the red part of the echelle spectrum of 
Procyon is insensitive to the iodine cell calibration and there- 
fore the search for stellar pulsations can be performed only by 
means of measurements of line equivalent widths (EW). The 
presence of telluric lines and limited number of strong lines 
sensitive to temperature (e.g. H a is located at the border of 
spectral orders and is therefore useless for accurate EW mea- 
surements) restricted our analysis to the spectral lines reported 
in Tabled 



4.1. Power spectrum 

Since the expected variations in the EW induced by the stel- 
lar pulsations are of the order of a few ppm, we adopted an 
accurate approach for the measurement of the EW, following 
Kjeldsen et al. I ( 1 1995b . Instead of fitting the line profile for de- 
termining the EW, we computed the flux in three artificial fil- 
ters, analogously to Stromgren Hg photometry. The first filter 
was centred on the considered line at the wavelength /In, with a 
band-width AAq, and the other two on the near continuum with 
respect to the line at Aq, one placed on the red region at A r , with 
band-width AA r , and the other on the blue one at Ab, with a 
band-width A/it,. We then computed the EW, or more precisely 
the line-index, in terms of the radiation flux measured in the 
selected filters T(AA): 



EW 



r(A^) 



A/in 



AA r + A/ib 



T(AA T ) + T(AA b ) 



(3) 



for two different central filter widths, AAi. and AA%, in order to 
obtain two values for the EW, EWi and EW2, from which the 
quantity EWb es t = EW'" +1 /EW"' can be constructed. The aim 
of this procedure is to construct a quantity that is less sensi- 
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time (JD) 



Fig. 13. Time series of EW measurements as derived from the 
combination of all the lines reported in Table [3] 



E S.Oxlo 4 

S 6.0x10 4 




Frequency [fiHzj 

Fig. 14. The weighted power spectrum EWb es t- < EW|, est > 
derived from the combination of the measurements of all the 
lines of Table [3] There is no excess of power in the range 0.5 - 
1 .5 mHz, but the background power in the range 0.2 - 0.7 mHz 
could be a signature of stellar granulation. 



tive to the continuum fluctuations, determining m so as to min- 
imize the scatter of the measurements in a similar way as de- 
scribed for the radial velocity analysis. The procedure was re- 
peated for all the lines reported in Table [3] and then, in order 
to improve the S/N, we combined all the EWbest's obtained for 
the examined lines, assuming for them the same sensitivity to 
temperature changes, and derived a weighted mean value. In 
Figs. [T3l and [T4l we show the time-series and the weighted PS 
of the EWbest, respectively. The weights adopted are w, = 1 /<x 2 , 
where the cr,'s are the standard deviations of the data. 




Frequency (mHz) 



Fig. 15. Smoothed power density spectra (PDS) of Procyon 
granulation noise as obtained from SARG (d ashed line, thi s 
proje ct) and AAT/ESO [EW 199 9. squares. dKieldsen et al 



Il999h l s pectroscopy and M OST dMatfhews et al.ll2004l) and 
WIRE dBruntt et al. I [2005) space photometry (continuous 
lines). On background the PDS derived from the PS of Fig. 
fT4lis shown. Power density is expressed in ppm 2 //Hz _1 . 

granulation. In this procedure we optimize the S/N in the PS 
by computing the weighted means of the measurements, by us- 
ing weights that minimize the mean noise for each line in the 
interval 0.2 - 0.7 mHz, and then rescale the weights for each 
line in such a way that the weights, wu satisfy the relationship 
Wi = 1 (3K, where J?l is the mean of the peak amplitudes in 
the range 0.2 - 0.7 mHz. 

In order to evaluate the power generated by the granula- 
tion of the star in the PS of the EW's and make a comparison 
with data obtained with different methods we need to define a 
quantity independent of the temporal length of the data, that is 
obtained by computing the power density spectrum (PDS), a 
measure of the power per frequency resolution element that is 
therefore independent of the length and sampling of the time- 
series: 



PDS(v) = A 2 (v) 



f 

Jo 



N(v)W 2 (v)dv 



-l 



(4) 



where A(v) is the amplitude (in ppm) in the PS, W(v) the win- 
dow function, and N(v) the Nyquist frequency. The result of our 
analysis is shown in Fig. [15] where the dashed line (this project) 
indicates the smoothed power density spectrum of Procyon 
granulation noise, and background the PDS of the PS shown 
in Fig. [H 



4.2. Granulation noise 

The combined power spectra do not show evidence of power 
excess in the range 0.5 - 1.5 mHz, as is clear from Fig. [14] 
However, at low frequencies, some concentra tion of power. 



simila r to the background power reported by Kjelds en et al. 



(1999) "or aCen A, an d attributed to s t ellar g ranulation, is ap- 
parent. Analogously to Kjel dsen et al. ( 19991) . here we attempt 
to give an estimate of the amount of power related to Procyon 



4.3. Comparison with other granulation noise 
determinations 

The literature reports three previous claims of g ranulation noise 



power detection in Procyon. The first is by iKjeldsen et al. 



(1999), who found a granulation noise slightly above the solar 
level in the frequency interval 1-2 mHz, by using the Balmer 
EW's. These determinations are shown in Fig. Q3] as squares 
(EW 1999). The other two claims, instead of spectroscopic, 
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concern space p hotometric measurem ents carried out by the 
satellites MOST dMatthews et al. 120041). operating i n the band- 
width 3700-7000 A. and WIRE feruntt et al. l2005t) . operating 
with the white light star tracker. 

While our measur ements are directly c omparable with the 
spectroscopic ones of IKjeldsen et al. I ( 1999b . in order to com- 
pare them with the photometric ones, we need to pass from 
intensities, /, to temperatures, T, and from the latter to EW's. 
This is accomplished by using appropriate sensitivity indices; 
for in tensity we used dlnl/dlnT = 4.5 (IKjeldsen & Bedding 



Table 4. Properties of the p-mode oscillations in Procyon. 



1995b . and for EW, dln(EW)/dlnT = -7.8, as deduced from 
the analysis of Bedding et al. (1996) for Fe I lines at T — 
6500 K. The results are shown in Fig. [15] where MOST and 



WIRE PDF curves are plotted together with IKjeldsen et al. 
(1999) data and our curve. 



It appears that, while WIRE and our distribution of power 
are mutually consistent, at least in the frequenc y interval 0.2 - 



2mHz, and corroborated by the few data of Kjelds en et al. 
(1999), VIOST shows a power density level higher by a fac 
tor of four or more. Uncertainties in sensitivity index calibra- 
tion of / and EW with respect to T may lead to some dis- 
crepancies, but not so large to explain such a significant dis- 
agreement between MOST and other data. In order to ren- 
der the MOST data comparable with ours, it should be | 
31n(EW)/(51nT | < 4.5 and dlnl/dlnT > 4.5. This is un- 
likely because the value of dln(EW))/<91nT adopted here has 
been derived by observations of th e same lines we use here at 
the same temperature of Procyon (Beddin g et al. 1996 ) and a 
valu e of d In I/dlnT > 4 .5 is in disagreement with the models 
of iBedding et al. I dl996b . This seems to imply that MOST did 
not detect granulation noise, probably because of some instru- 



menta l or stellar spurious signal as suggested by Bru ntt et al. 
d2005l) . 

Differently f rom th e first spectroscopic determinations of 



IKjeldsen et al I (119991) relative to a few points in the high 
frequency side of the granulation noise PDS, our EW deter- 
minations constitute the first spectroscopic measurements of 
Procyon granulation in the frequency range 0.1 - 2.0 mHz, so 
providing an independent measure of and establishing an upper 
limit to granulation noise power. The fact that at low frequen- 
cies, below 0.2 mHz, our power level is higher than that shown 
by WIRE might be caused partially because we did not correct 
the WIRE data at low frequencies for the effect of high-pass fil- 
ter and partially because of the possible effect of stellar activity 
in Procyon that affects spectral lines but not the continuum. 

5. Conclusions 

The analysis of the PS of Procyon obtained by SARG high 
resolution spectrograph shows unequivocally the presence of 
a solar-like p-mode spectrum as clearly demonstrated in Fig. 
[7] The main seismic properties of Procyon are summarized in 
Table E] 

We identified with a good margin of certainty 1 1 individual 
p-mode frequencies with / = 0,1,2 and 7 < n < 22 in the 
frequency range 500 - 1400 //Hz, as reported in TableQ] which 



Large separation Av 


55.90 ±0.08 //Hz 


Small separation 6vq2 


7.1 ± 1.3 //Hz 


Asymptotic relationship e 


1.913 ±0.025 


Frequency scatter per mode 


1.23 //Hz 


v(Z = 0, n = 16) 


1001.3 ±0.4//Hz 


Mean amplitude per mode (/ = 0, 1) 


0.45 ± 0.07 ms- 1 


Mode lifetime 


2.0 ± 0.4 d 



The large frequency separation deduced from the fit of 10, 
out of the 11 identified frequencies, to the asymptotic rela- 
tionship (HJ is consistent with the results from the comb re- 
sponse analysis, both individual and cumu lative, and its value 



agrees fairly well with that determined by Eggenberg er et al. 



( 2004), Av = 5 5.5 //Hz, but shows s ome discrepancy with that 
determined bv iMartic et al. I d2004 Av = 53.6 //Hz, though 
in a previous article oMartic et akl IT999) the authors found 
Av = 55.6 //Hz, as is evident from their Fig. 12. 

The large frequency separation we have determined is con- 
sistent with an evolutionary model of Procyon, constructed 
with solar chemical abundance, elemental diffusion and mass 
loss, having a mass of 1.48 M Q , an age of 1.8 Gy, a lumi- 
nosity of 6.9 L , a radius of 2.04 R , and a surface convec - 
tion z one deep the 8 % of the stellar radius dBonanno et al. 
20061) . It is in very good agreement with other theoretical stud 



Eggenberger et al. 



ies of Procyon A (Barba n et al. 1119991; IChabover et al. 111999 



2005) 



are lar gely consistent with those found by Eggenberg er et al. 
(2004)? 



The small frequency separation we determined from one 
I - 2 and t hree I — \ modes is larger by 1 .6 <x than that de- 
termined by Marti c et al. ( 2004h , 6vq2 — 5.1 //Hz, and that we 
derived from Eggenberger et al. I (120041) data, <5vo2 = 4.9 //Hz, 
but our value is not strongly constrained by our measurements, 
as only a few modes with / + were identified. 

The mean amplitude per mode, for modes with / = 0, 1, is 
about 0.45, twice larger than the solar one. 

The mode lifetime has been determined by means of two in- 
dependent methods whose results are mutually consistent and 
give a lifetime of about 2 d. This short lifetime rules out the 
possibility of a coherent, over-stable excitation of modes in 
Procyon, but indicates that oscillations are excited stochasti- 
cally by turbulent convection as in the case of the Sun. 

We used the red part of the observed spectrum for deter- 
mining, through the measurement of the EW's, an upper limit 
to granulation noise pow er that resulted to be in agreement with 
previous spectroscopic ([Kjeldsen et al. II 19991) and space pho- 
tometric ( iBruntt et al. 112005b measurements, but in significant 
contrast with MOST data (Matthe ws et al. Il2004 . 

Our radial velocity data are available upon request to the 
first author of the present article. 
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